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Abstract 

We construct a spin foam model of Yang-Mills theory coupled to gravity by 
using a discretized path integral of the BF theory with polynomial interac- 
tions and the Barret-Crane ansatz. In the Euclidian gravity case we obtain 
a vertex amplitude which is determined by a vertex operator acting on a 
simple spin network function. The Euclidian gravity results can be straight- 
forwardly extended to the Lorentzian case, so that we propose a Lorentzian 
spin foam model of Yang-Mills theory coupled to gravity. 
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1 Introduction 



Spin foam (SF) models of quantum gravity represent a promising approach 
for formulating a consistent theory of quantum gravity [1, 2]. It is a dis- 
cretized spacetime path-integral approach, and the main result is the finiten- 
ness of the 3-geometry to 3-geometry transition amplitude for a non-degenerate 
space-time triangulation [3]. Incorporating matter in the SF framework has 
been started in [4]. However, the constructions given there were algebraic, 
without a direct connection to the discrete path-integral considerations. The 
difficulty in the path-integral approach is that the usual matter fields couple 
to gravity via the fierbein one-forms, while in the SF framework the funda- 
mental field is a 2-form B which can be considered as an exterior product 
of two fierbein one-forms. Still, in some cases the coupling of matter can be 
expressed conveniently in terms of the B fields, like in the Yang-Mills (YM) 
theory case [4] 



where Fki = e%e"F^ u and is the inverse of the fierbein matrix e a . F is the 
YM two-form field strength, and Tr is the gauge group trace. 

In such cases one can use the discretized path-integral techniques devel- 
oped for the SU(2) BF theory by Freidel and Krasnov [5], in order to derive 
the SF amplitude. Note that Oritti and Pfeiffer have proposed recently a 
YM euclidian SF amplitude without using the path- integral approach [6]. In 
order to better understand their result, we will study the discretized path- 
integral for a general BF theory with polynomial in B interactions by using 
the Friedel-Krasnov approach. The algebraic structure of the correspond- 
ing SF amplitudes is such that it is straightforward to adapt them via the 
Barrett-Crane ansatz [7, 8] to the case when the B field is a simple bivector, 
i.e. when it describes a realistic spacetime geometry. 

In section two we generalize the Friedel-Krasnov results to the case of 
an arbitrary Lie group G. In section three we perform the Barrett-Crane 
reduction of the results from the section two. In section four we construct 
the YM spin foam amplitudes. In section five we present our conclusions. 




(1) 
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2 Path-integral approach 



Consider a BF theory with a polynomial interaction, depending only on the 
B field 

S = [ {B A F) + [ P(B) d A x , (2) 
Jm Jm 

where M is a four- dimensional manifold and B is a two-form taking values 
in the Lie algebra g of the Lie group G. F = (gL4 7 + ff K Aj A A^T 1 is the 
curvature two- form taking values in g and (, ) is an invariant quadratic form 
on g. 

We would like to define the path integral 

Z = J VAVBe iS , (3) 

via the discretization procedure provided by a simplical decomposition of M. 
We will use the generating functional technique, so that we introduce 

Z[J]= J VAVBe iS+i ^ BAj ) (4) 

and 

Z [J]= J VAVBe^M {BA{F+J)) . (5) 

Z[J]=exp(,/ M P(-4))z [J] , (6) 
and Z = Z[0]. 

In order to define Z let C be a simplical complex associated to M, and let 
C* be the dual complex. We denote by a, r and A a 4-simplex, a tetrahedron 
and a triangle of C respectively, and by v, I and / the corresponding dual 
cells, i.e. a vertex, a link and a face, so that 

v = a* , l = r* , / = A* . (7) 

Let B field be piece-wise constant in the 4-polytopes formed by the pairs 
(A, A*) such that the only non-zero components of B lay in the A planes 1 ". 
Then 

S = J2(B A ,F f )+ J2 C I ... I ,(A,...,A')B I A .-.B I A , , (8) 



Hence 



A,..., A' 



TThis discretization procedure is more suitable for treating interactions than the one 
used in [5] , which was based on the delta-function B fields. See [9] for basic notions about 
the forms on simplical complexes. 
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where 

Ba= L B ■ F '=l F ■ (9) 

Given the variables (9), we define 

Z [J]= fU dA iU dB Ae l ^ {BA ' F,+Jf) , (10) 

J l A 

where A[ = j l A. Following the ref. [5] we will fix the measures of integration 
dA and dB such that 

w = /n^n*fnM" ) )) » ( n ) 

J i f \ief v 7 J 

where g x = e Al , and are the Lie algebra elements associated to the vertices 
v of the face / such that 

YlJ? = Jf ■ (12) 

vef 

The formula (11) differs from the one given in [5] by the absence of the factor 
Ilvef P(Jf V ^)- We do not put this factor because the simpler expression (11) 
is also invariant under the gauge transformations 

9i^>h v gih~} , J f (v) -> h v J f (v) h' 1 , h v ,h v >eG, (13) 

where v and v' are the ends of the link I [5] . Also when J — > the expression 
(11) gives the standard formula for the partition function for the BF theory 
[10, 1]. Hence we take the expression (11) as the definition of Z (J). 
The expression (11) can be further simplified by taking 

jM = J^) = ... = ^ , (14) 
n f 

where nf is the number of vertices of the face /. By rescaling Jf — > rifJf we 
can get rid of the nj factors. By using the group theory formulas 

6(g) =]TdimAxA(3) , (15) 

A 

and 

f da D (Al) (n\ D (A2) (n) D iA:i} (a) D (A4) (a) - V /^i^sMO / r yA 1 A 2 A 3 A 4 W\ * 

(16) 
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where is the group representation matrix in a representation A, and xa 
is the corresponding trace, as well as the associated graphical calculus [5, 1], 
one arrives at 

2o=eii dimA / nca,)n*6(A / ( W ),t I(1 , ) , e j /w) , (17) 

A f M f I v 

where ^5 is the spin network function associated to the pentagram graph 
and C(i) is a constant coming out from the recoupling theory calculus. 

The problem with (17) is that the sum over the irreps diverges. This 
sum can be regularized in the J = case by passing to the quantum group 
G q = U q (g), q root of unity [11]. One could generalize this procedure for 
the non-zero J case, but since in this paper we are not interesting in the 
topological gravity case, we will not elaborate upon this. 

3 The Barrett-Crane ansatz 

General relativity can be understood as a BF theory with constraints since 

J y/te\Rd*x = J e abcd e a Ae b AR cd = J B ab A R ab , (18) 

where e a are the fierbein one-forms and R ab is the spin-connection curvature 
two-form. Hence G = SO (4) in the Euclidian case or G = SO (3, 1) in the 
Lorentzian case and the constraint is B ab = *(e a A e b ). In the SF framework 
this constraint is implemented as a restriction on the set of irreps {A/} such 
that 

(T(A / ), *T(A / )> = e abcd T ab (A f )T cd (A f ) = , (19) 

where T(A) are the generators of G in the representation A [7, 15, 8]. The cor- 
responding solutions are called simple irreps, which we denote as N. Equiv- 
alently, N is the irrep containing an 5*0(3) invariant vector [13]. 

Given the set of simple irreps {Nf} one postulates that Z is of the same 
form as in the topological case 

Z = Y.UMNf)ilA 1 (N m )HA (N m ) , (20) 

N f f I v 

but now the amplitudes Ai are different. In the Euclidian case = 

j G |Z, and the face amplitude A 2 = dim N = (2j + l) 2 [7]. In the Lorentzian 
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case A 2 = p 2 dp, where N = (0,p), p > [8]. The vertex amplitude A , which 
is associated to the pentagram spin net, is given now as 



. 5 10 

A (N ir --,N w ) = / Hdx v l[K Nl (x v (i),x v '(i)) , (21) 

JG / H v=l 1 = 1 



where x v are the points in the homogeneous space G/H and H = SO(3) 
[7, 12, 13, 8]. The propagators K are given as 

K N (x,y) = (0\D^(g x g; 1 )\0) , (22) 

where |0) is the invariant vector from N, which corresponds to the identity 
irrep of the subgroup H contained in N . There is a freedom in choosing the 
edge amplitude A ± , and a choice which makes the sum (20) finite is the spin 
net amplitude for the theta-four graph [14, 3]. 

Given these results and the formula (17) for the topological case, it is 
natural to propose the following expression for Z (J) in the non-topological 
case 

z °(j) =En^( iv /)n^( iv /(o)n $ 5(iv / (,);e^)), (23) 

N f f I v 

where 



.5 10 

<5> 5 (N 1 ,---,N 10 ;e J \---,e J ™)= / ]J dx v ]J K Nl (x v{l) , x v , {l) - e Jl ) , (24) 



5 10 
,G / H v =l l=i 

and the propagator 

K N (x,y;e J ) = (0\D( N \g x e J g?)\0) , (25) 

is the "source" propagator introduced in [13]. Hence <3> 5 is an example of the 
simple spin network functions introduced in [13]. 



4 The Yang-Mills model 

The YM action on C is given by [16] 

A(A 



\ 2 S YM = E-^(ReTrU(A)-n) 



(26) 
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where U is the triangle holonomy associated to the gauge group G and n is 
the dimension of the matrix U . A(A) is the area of the triangle A, while 
A(A*) is the area of its dual face. Let us rewrite this action as 

\ 2 S YM = E 6V £ff (Re Tr U(A) - n) , (27) 

where V(A, A*) = |A(A)A(A*) is the 4-volume of the 4-polytope (A, A*) 
[9]. This can be further rewritten as 

^ym ^EE 6 ^ a ^ (ReTrU(A)-n) , (28) 

which is suitable for the spin foam formalism. 

The Oriti-Pfeiffer proposal [6] can be understood as an approximation 

V(A,A*)« £ C(A,A'){B A ,*B A ,) . (29) 

This is an approximation because V(A, A*) depends in general on the B 
fields from the adjoint 4-simplices which share the polytope (A, A*), so that 
for an exact formula one would have to include the terms with A' e o' . A 
further approximation is to take the C's to be constant up to the orienta- 
tion sign factors sign(A, A'). This happens for the symmetric lattices when 
V(A, A*) = (2/5)V(a), where V(a) is the 4-volume of the 4-simplex a. Since 

V W = QfTTi E sign(A, A')(5a,*5a') , (30) 

6U ' 4 ' A, A' 

one obtains 

C(A,A')«^sign(A,A') . (31) 
One can then define 

z YM = / n d 9e n d 9i n dB ± e^-^™ = / n ^ e , (32) 

where e's are the edges of the complex C. By using the previous results we 
obtain 

zym « e n ^ ( jv/) n ^ (%o ) n $ 5 (%, } ; e >i (33) 

Nf f I v 



6 



where (3 is a numerical constant and 

(ReTrU(A) - n) I Q q \ 

*m j ^ (A ' A,) W'W ■ (34) 

A natural further simplification is to take 

A(A) = y/j{j + 1) , (35) 

in the Euclidian case [6]. In the Lorentzian case the labels of the triangles 
can be related to the areas of space-like triangles [8], and hence it is natural 
to take 

A(A) = V? + I • (36) 

One can show that 

e-^ 5 (N f{v) ; e '/<->)] = £ C£::£ ( ai • • • a^e'^^ ■ ■ ■ f3 10 ) , 

(37) 

where S(T) is the operator obtained from (34) by replacing the <9J's with 
the T(jV)'s. It acts in the space = (8>ago- V(Nf)- The coeficients C are 
given as products of five Barrett-Crane intertwiners 

«=EC^ W • (38) 

N 

From (37) it is clear that in the Euclidian case the expression e~ lf3Sv $5 1 j = o 
is well defined since TV's are finite-dimensional and hence all the sums are 
finite. The expectation value in (37) is the matrix element of an exponential 
of a finite matrix, which exists. 

In the Lorentzian case, the operator representation (37) is not that useful 
for showing that e~ ll3Sv ^^\ j=q is well-defined. We expect it to be well defined 
because in the Euclidian case the propagators* 

1 x sin(2? + l)d(x, y) ,, s 

Kj(x,y)= } J / ; ,y \ , cosd(x,y)=x-y , 39 
(2j + 1) smd(x, y) 

get deformed when the sorce J is turned on as 

cosrf — > cosrf cos I J\ + sinjrf sin | J| , (40) 

*The extra factor (2j + follows from the formula (22), since the SU{2) invariant 
vector has a normalisation factor (2j + l) -1 / 2 . 
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where 

sin j d = cos 9 sin cos 0j + sin 9 cos cos 9j (41) 

for a; = (cos 9, sin 0, 0) and y = (cos 0, sin cos a, sin sin a, 0). The angles 
9 j and 0j are determined by the 3- vectors x and J and y and J, respectively, 
which are associated to the corresponding SU(2) subgroup elements. The 
deformation (40) is qualitatively similar to the deformation d — > d + \ J|, so 
that it is no surprise that the integral (24) is finite, given that the J = 
integral is finite. Since the Lorentzian propagator 

T . , s sinpd(x,y) , . . in . 

K ^ y) = psinhd(x,y) ' ™ hd ^y)=*-y > ( 42 ) 

can be obtained from the Euclidian one via analytical continuation 2j '• + 1 — > 
ip and d — > id, we then expect to have a deformation 

cosh d — > cosh d cos | J| + sinhj d sin | J\ , (43) 

where 

sinhj d = cosh # sinh cosh 0j + sinh 9 cosh cosh #j (44) 

for x = (cosh 9, sinh 0, 0, 0) and y = (cosh 0, sinh cos a, sinh sin a, 0). This 
is similar to the deformation d — > d + i|J|. Since the J = integral (24) 
converges, we then expect that the J ^ integral will converge as well. 



5 Conclusions 

The discretized path-integral approach to spin foam amplitudes based on the 
generating functional technique is very useful, because the algebraic structure 
of the expressions one obtains in the topological case is such that it can 
be straightforwardly extended to the non-topological cases. This algebraic 
structure is the tensor category of irreducible representations of the symmetry 
group, so that the spin foam transition amplitudes can be understood as 
functors constructed from the special morphisms, which are the spin networks 
[17]. Even in the non-topological case, when the simple irreps iV do not form 
a tensor category, one is actually dealing with the tensor category of the 
irreps of the subgroup H [17], which in the YM case are just the trivial 
identity irreps, so that it is natural that we have proposed the generating 
functional in terms of the simple spin network functions. 
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Because of this, we consider our proposal for the vertex amplitude (37) 
more appropriate for the spin foam formalism than the one proposed in [6], 
which is 

H f „ ) - = J2(^H v r 1 (a 1 ---a w \e-^ s ^\a 1 ---a w ) . (45) 

dimii, — * 



V 



Although this is a reasonable proposal, this expression does not have a simple 
interpretation in terms of spin network evaluations. Consequently, it is not 
clear what would be a Lorentzian extension of this expression. On the other 
hand, our expression (37) is based on a simple spin network evaluation, and 
hence it makes sense in the Lorentzian case. 

We expect that the YM amplitude (33) is finite both in the Euclidian 
and in the Lorentzian case because of the very rapid convergence properties 
of the J = sums in the Euclidian case [18]. 

Note that the structure of the vertex amplitude is such that one can use 
the perturbative expansion for small (3 

e -^ $ i = r (1 _ ij3 § v + . . . )$ i = Ao _ i(3 r^i , . . _ (46) 



J=0 



J=0 



J=0 



The second term in (46) can be interpreted as the pentagram spin network 
amplitude with the insertions of two endomorphisms T(N) and T(N') at 
the edges N and N'. If we neglect the higher-order terms in the expansion 
(46) we will obtain a spin foam amplitude which is a sum of amplitudes 
where a string of pentagram spin networks is connected by the matter edges 
carrying the adjoint representation of G. This amplitude is of the general 
type proposed in [4], but it is more complicated then the ones considered 
there. In [4], only the amplitudes where the matter edges link the vertices 
of the pentagrams were considered, while the amplitude corresponding to 
(46) has the matter edges linking the edges of the pentagrams and there are 
T(N) insertions. This demonstrates that in order to understand completely 
the coupling of matter in the spin foam formalism, it is necessary to study 
the corresponding discretized path integral. 
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